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Abstract 

Stochastic phenomena in which the noise amplitude is proportional to the 
fluctuating variable itself, usually called multiplicative noise, appear ubiqui- 
tously in physics, biology, economy and social sciences. The properties of 
spatially extended systems with this type of stochasticity, paying special at- 
tention to the non-equilibrium phase transitions these systems may exhibit, 
are reviewed here. In particular we study and classify the possible univer- 
sality classes of such transitions, and discuss some specific physical real- 
izations including depinning transitions of non-equilibrium interfaces, non- 
equilibrium wetting phenomena, synchronization of spatially extended sys- 
tems and many others. 

1 Introduction: Multiplicative Noise 

Since the pioneering description of Brownian motion as a stochastic process, Lan- 
gevin representations of statistical systems have played a central role in the devel- 
opment and understanding of statistical mechanics Q|2][3]]. The introduction of 
noise on otherwise deterministic (noiseless) systems have proven to be essential 
in order to properly describe numberless problems in physics, chemistry, biology, 
economy or social sciences. In some cases noise may lead to rather surprising and 
counterintuitive behaviors. A well known example of this is stochastic resonance 
in which the response of a stochastic system to an external periodically oscillating 
field is enhanced by increasing its internal noise amplitude [4]. Other instances 
are resonant activation [5], noise-induced spatial patterns H3|7), and noise induced 
ordering transitions (NIOT) (8l|9| to name but a few. 

Most of the first tackled examples, including Brownian motion, were prop- 
erly described by Langevin equations with additive noise, i. e. state-independent 
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fluctuation amplitudes. However, in order to properly account for more complex 
situations multiplicative noise (MN), i. e. a variable-dependent noise amplitude, 
might be (and, actually, is) required. A prototypical example of this are stock mar- 
ket fluctuations in which the dayly variation of a given company quotation is given 
by a certain percentage of its initial value, i. e. variations are larger for larger start- 
ing prices. Another example are time fluctuations in the number of molecules of a 
given species present in a catalytic reaction: the fluctuation amplitude depends ob- 
viously on the number of molecules at each time. In the same way, fluctuations in 
population dynamics models as the Lokta-Volterra can be modeled by multiplica- 
tive noise [10]. Also, situations in which some control parameter fluctuates may 
generically lead to the presence of MN l6l ll II . Some more example are the equa- 
tions describing the fluctuations of advancing fronts fT2ll . or Schrodinger equations 
in the presence of disordered potentials [ 13 1. 

It is well established by now that many of the situations in which noise gener- 
ates counter-intuitive features can be modeled by Langevin equations with MN. As 
a rule of thumb, in systems with additive noise, the smaller the noise amplitude, the 
more "ordered" the stationary state. On the other hand, for systems with MN this 
is not necessarily the case, and, as said before, a wealth of instances have been de- 
scribed where the degree of ordering is enhanced upon enlarging noise amplitude; 
noise induced phase transitions are a good example l8ll6ll9l [T4l[T5lfT6l . 

Before proceeding further, let us stress that among many other applications, 
Langevin equations (with or without MN) have turned out to be a highly conve- 
nient theoretical tool for the analysis and characterization of critical phenomena, 
including both static and dynamical aspects; an application that will be of special 
interest along this paper. For instance, models A and B (both of them defined by 
Langevin equations [3 1) have become paradigmatical examples describing the uni- 
versal features of Ising-like phase transitions, with non-conserved and conserved 
dynamics respectively, in the simplest possible way 1 . 

In this review we focus on phase transitions in spatially extended systems with 
multiplicative noise. We will explore their phenomenology, paying special atten- 
tion to their critical properties, and will review some applications to physical sys- 
tems, including pinning-depinning, wetting problems and synchronization transi- 
tions, among many others. 

1 Other symmetries and/or conservation laws can be easily implemented in this language, gener- 
ating a zoo o Langevin models: C, D,... H. |3J 
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2 Multiplicative Noise basic equation 



In has been recently shown [ 16 1 that the following multiplicative noise (MN) Lan- 
gevin equation intended in the Stratonovich sense 2 is the minimal model for 
noise-induced ordering transitions in spatially extended systems iffl rffl : 

d t (f>(x, t) = -a(j>(x, t) - bcjf^x, t) + DV 2 4>(x, t) + a(j)(x, t)rj(x, t). (1) 

rj(x, t) is a Gaussian white noise, 4>(x, t) the order-parameter field that usually 
represents a density or an activity-degree field, a, b > 0, D > and a > are 
parameters, and p is an integer constant obeying p > 2 (typically p = 2 or 3). As 
discussed in [ 16 1, for this equation the larger the noise amplitude the more ordered 
the stationary state. A very similar equation has been proposed by Pikovsky and 
Kurths as a simple stochastic description of synchronization transitions (ST) in 
extended systems |H71[TH1 3 (see also section [SJl. 

Eq.Q exhibits, among other non-trivial features, a non-equilibrium phase tran- 
sition as the control parameter a is varied, that has been profusely studied lfT9ll20l 

E3- In particular, for a > a c , where a c corresponds to a transition point, the 
only stationary state is 4>(x) = 0, defining an absorbing phase 4 . On the other 
hand for a < a c the system seats in an active phase with (<p) ^ 0. This defines an 
absorbing-state type of phase transition. 

Even if, strictly speaking, a MN Langevin equation is any in which the noise 
is a certain function of the stochastic variable, in what follows, we restrict the ad- 
jective multiplicative to describe a noise linear in the field, as the one in Eq.Q. 
This excludes explicitely, among many other cases, the Reggeon field theory Lan- 
gevin equation in which the noise amplitude is proportional to the square-root of 
the variable (see below). 

Let us start the analysis of Eq.Q, the most basic MN equation, 5 by describing 
its three contributions: 

1. A single-site deterministic force, F(x, t) = —acj){x, i)—b<j?(x, t), which can 
be derived from the potential V[(f>] = |(/> 2 + ^i<p p+1 - This is the simplest 
potential describing a continuous bifurcation (transition) from a <fi = state 

2 See Appendix A for a brief discussion of the Ito versus Stratonovich dilemma. 

3 To the best of our knowledge this is the first time that Eq.Q appeared in the literature. 

4 Sometimes we will refer to this as "quasi-absorbing" phase for reasons that hopefully will be- 
come clear in what follows. 

5 The effect of MN has been studied also in other type of Langevin equations, as for example 
the Swift-Hohenberg one |7||23I . that has become in paradigm in the study of pattern formation, or 
equations with conserved order-parameter as the Cahn Hilliard one 1 23 6 1. We will not explore these 
cases in the present review. 
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(for a > 0) to a non trivial state, cf> ^ (for a < 0) 6 . Higher order terms 
could be included in V(h) but, as it is well known in the folklore of critical 
phenomena, they are expected to be irrelevant in what respects universal 
properties, at least as long as b remains positive, ensuring stability. 

If for symmetry reasons the order parameter of the problem under consider- 
ation exhibits a Z2 inversion (or plus/minus or up/down) symmetry then an 
odd value of p (typically p = 3) is required, 7 while p = 2 will be used by 
default in the absence of such a symmetry as the most basic choice. 

2. A diffusion term, DV 2 <ft(x, t), which, as usual, is the force derived from a 
surface tension D J d d x^(V<p) 2 tending to smooth away fluctuations. Higher 
order terms in a gradient expansion could also be introduced without altering 
the critical properties. 

3. A noise term proportional to <fi(x, t), i. e. a MN. 

In order to put the problem under perspective, and for pedagogical reasons, let 
us compare this Langevin equation to other very well-known ones: 

(A) Model A describing Ising-like phase transitions with non-conserved dy- 
namics EElEa: 

d t <j>{x, t) = -aqb- b(p 3 + D S/ 2 cp(x, t) +arj(x,t). (2) 

Apart from the cubic term imposing up-down symmetry, the only difference with 
Eq.0 is in the noise amplitude, constant in this case. 

(B) The Reggeon-Field theory (RFT) known to describe phase transitions in 
a vast class of systems with absorbing states (directed percolation (DP) and the 
contact process being paradigmatic instances l26ll27ll28ll29ll30ll3Tl ). 

d t cp(x, t) = -a(p - bcj) 2 + DV 2 (p(x, t) + cr^4>(x, t)r)(x, t). (3) 

The only difference with Eq.Q is in the noise amplitude, proportional here to the 
square-root of the activity field. The difference between the absorbing states in 
this class and the "quasi-absorbing" states of the MN family is discussed in detail 
in [32] (see also Appendix B). In a nutshell the difference is that "true" (RFT) 
absorbing states correspond to a non-integrable singularity at the origin (<ft = 0) of 
the single-site probability distribution function (present even in the active phase) 

6 This is usually called a transcritical bifurcation in the literature of dynamical systems. 
7 Observe that if p = 3 the order parameter can be either positive definite, or defined in the whole 
real axis. This should not affect critical properties. 
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while such a singularity does not exist or is integrable in the active phase of 
problems in the MN family. 

The three equations presented so far, with noise amplitudes <p°, and (j) 1 
respectively, represent three main groups of phase transitions occurring in systems 
without anisotropics, long-range interactions, extra components, non-Markovian 
effects, or quenched disorder. An ever-growing number of critical phenomena in 
very different contexts have been reported to belong to any of these very robust 
three classes. While for both the Ising class (Model A), and the directed percolation 
(RFT) classes several comprehensive reviews exist in the literature (see for example 
0J|23|28j|29|[31 1) to the best of our knowledge there is no compilation of known 
results for the MN equation and variations of it. It is the purpose of this review 
to report on the interesting and rich phenomenology of this and closely related 
systems, as well as on their many physical realizations. 

3 Connections with related problems 

Before analyzing the properties of Eq.Q, for the sake of generality, we enumer- 
ate now its connections with other relevant problems in non-equilibrium statistical 
physics. In particular it is related to: 

1. The Kardar-Parisi-Zhang (KPZ) equation 

The KPZ equation describes genetically the roughening of interfaces under 
non-equilibrium conditions and has attracted an enormous deal of attention in the 
last decade I33ll34ll35ll . To see how the connection between MN and KPZ works, 
Its connection with the MN equation comes out by performing the so-called Cole- 
Hopf transformation 4>(x, t) = exp[h(x, t)] to Eq. 

d t h(x, t) = -a-b e (p_1) h + DV 2 h + D{Vhf + ar](x, t) (4) 

which is the celebrated KPZ equation, describing an interface of height h{x,t) 
(with respect to a reference substrate at h = 0) in the presence of an additional 
bounding term, —b e^ -1 )^, pushing regions where h(x,t) > towards negative 
values. The first two terms can be written as derived from a potential V(h) = 
ah + ^ b ^ e[( p ~ 1 ) fc ] (observe that a controls the potential slope for large values of 
\h\). Eq.© is usually called in the literature a KPZ with an upper wall. If instead 
we had defined 4> = exp(— h), the sign of the emerging KPZ-non-linearity would 
be negative, inducing a lower wall (i. e. the argument of the exponential bounding 
term changes sign, and therefore negative values of the height variable are pushed 
towards positive ones). In other words, a KPZ interface with positive non-linearity 

8 This means that the only true stationary solution is a Dirac-delta distribution at <f> = 0. 
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in the presence of an upper wall is completely equivalent to a KPZ with negative 
non-linearity and a lower wall. These two cases can be mapped one into the other 
and are described by Eq.Q. We will refer to this as Multiplicative Noise 1 (MN1) 
class. 

The active phase, (0) / 0, is translated in the interface language into a pinned- 
to-the-wall phase while the absorbing one maps into a depinned phase (with in- 
terfaces escaping towards infinity or minus infinity at a constant average velocity). 
The absorbing-active phase transition becomes a pinning-depinning one within this 
language. 

The curious reader could pose himself/herself the following question: what 
happens if we consider a KPZ interface in the two remaining cases, i. e. with pos- 
itive non-linearity in the presence of a lower wall, or with a negative non-linearity 
and an upper wall? Can these problems be described in terms of a MN equation? 
Taking the following KPZ-like equation with an upper wall ( b e^ p_1 ^) and a neg- 
ative KPZ-nonlinearity, 

d t h(x, t) = -a-b e ip ~ 1)h + DV 2 h - D{Vhf + arj(x, t) (5) 

as starting point, and performing the inverse Cole-Hopf transformation we arrive 
at: 

d t Mx, t) = -acj) - b<j? + DV 2 cf)(x, t) - 2D ^^ X,t ^ + a<p(x, t)ri(x, t) (6) 

(p(x,t) 

which is identical to Eq.Q except for an extra singular term. Eq.@ is a modified 
MN equation, that we will name Multiplicative Noise 2 (MN2) equation. Obvi- 
ously the same result is obtained for a lower wall and a positive KPZ non-linearity 
by using <fi = exp(-h). As we will illustrate in the forthcoming sections the extra 
term —2D ^fj- x 'j)^ , which can also be written as 

- 2D(V(p(x, t)).(V \n(<p(x, t))) = ±2D(V(p(x, t)) ■ (Vh(x, t)) (7) 

(apart from being singular) introduces a relevant perturbation making MN2 essen- 
tially different from MN1 9 . 

Summing up: for a fixed value of the KPZ non-linearity the introduction of 
limiting upper or lower walls lead to essentially different physics. For a fixed wall, 
changing the sign of the non-linearity also changes essentially the behavior. 

Let us illustrate the physical reasons behind this fact. For that, we fix the sign 
of the non-linearity to be positive. For an upper (lower) wall the non-linearity pulls 

'Observe that the presence of Vh(x, t) in the equation for <f>'s Eq.0 means that MN2 finds a 
more natural representation in terms on h's than in terms of </>'s. 
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Figure 1 : Schematic representation of KPZ interfaces with a positive non-linearity 
in the presence of upper (MN1) and lower (MN2) walls respectively. The walls 
are depicted as rigid lines (corresponding to p — > oo; for finite values of p the 
walls become exponential cutoffs). Observe the difference between the two cases: 
while in one of them the constant driving force pulls the interface away from the 
bounding wall (MN1), in the other, the interface is pushed to the wall (MN2). 



(pushes) the interface from (against) the wall, while the constant term, acts in the 
opposite direction in order to keep the interface pinned. These are the two only 
mechanisms driving the interface average velocity, which imposing stationarity is 
given by 

(8 t h) = =±a±(b e Th ) + (D{Vh) 2 ). (8) 

At the critical point (e Th ) = and, therefore, a has to compensate (D(Vh) 2 }. Let 
us now consider a flat interface ((V/i) 2 = 0) at the critical point. While in one 
case the constant drift, a, pulls it from the wall, causing an effective repulsion from 
it, in the other one the flat piece is pushed against the wall, being therefore stable 
(see figure This mechanism leads to essentially different stationary probability 
distributions in both cases, and ultimately to different universality classes [21 1. 

The limit p — > oo describes a hard (impenetrable) wall: the interface cannot 
cross the h = limit in either the upper (h(x,t) < 0) nor the lower wall case 
(h(x, t) > 0). As we will show the "hardness" (or rigidity) of the wall is not 
a relevant feature. Indeed, p can be reabsorbed from the KPZ-like equation by 
modifying the Cole-Hopf transformation to a more general form: <fi = exp(m/i), 
where m is a conveniently chosen constant. 

2. Directed polymers in random media (DPRM) [33 1. 
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The connection with directed polymers in random media has also been known 
for some time. Indeed, the partition function of the DPRM can be easily seen to 
obey the following equation (imaginary time Schrodinger equation) [ 33 1 

d t Z(x, t) = [DV 2 + r](x, t)]Z(x, t) (9) 

including both: an elastic term (diffusion) and a random potential. Therefore any 
problem in which the possible values of the DPRM free energy, F oc In Z are 
restricted, owing to some unspecified reason, from above or from below, can be 
mapped into a MN equation, or into the KPZ equation with an upper or lower 
wall respectively. In particular, it has been pointed out by Hwa and coauthors that 
algorithms used routinely in molecular biology for sequence local-alignment are 
straightforward realizations of DPRM with a free-energy constraint and are related 
to the KPZ-problem with a lower wall [36||^. Other problems in condensed mat- 
ter physics, in which a Schrodinger equation appear including a random potential, 
can be ascribed to this class; this can include problems related to Kerr dielectric 
guides, propagation of optical pulses, multi-channel optical fiber networks, etc. 
See [ 13 1 and references therein. 

3. Noisy Burgers equation of turbulence ll37ll33l . 

Considering a KPZ equation and transforming v(x, t) = S/h(x, t) one obtains 
the (noisy) Burgers equation describing a randomly stirred vorticity-free (i. e. V x 
v = 0) fluid 

^ = d t v + v ■ Vv = 2V 2 v + Vri (10) 
at 

where we have fixed D = 2 to ensure invariance of the total derivative under 
rescaling [33 34 1. The presence of an extra upper wall generates an additional 
term which cancels out high velocities. This velocity constraint can mimic the role 
of viscosity in some problems in fluid-dynamics. 

4 Properties of Multiplicative Noise. The MN1 class 

In this section we explore the phenomenology associated with MN1 by employing 
different levels of approximation and tools. 

4.1 Zero dimensions 

An exact solution of the one-variable (zero-dimensional) case was provided by 
Graham and Schenzle [38] (see also [8| and [30 1 for a review). Indeed, the sta- 
tionary solution of the associated Fokker Planck equation can be easily worked-out 
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(see Appendix B) 



Let us just underline two curious non-trivial properties of this solution [38 1: 



-2bc/)i 



(11) 



• All the moments scale in the same way(!) For example, for p = 2 it is just a 
matter of algebra to obtain 



2m/2p 



m + a/a 2 



-a/2a< 



(12) 



m ~ |a— a c | /3m is equal to unity 1) for 



and, therefore, j3 m defined by 
an arbitrary p) for all m [38] El- This is a rather extreme case of anomalous 
scaling, having its roots in the fact that the asymptotic behavior is controlled 
by large rare events I39l l40l . 



The susceptibility (response to an external field) diverges not only at the crit- 
ical point, but in a finite interval in the absorbing phase, and the associated 
critical exponent varies continuously 1191 . 



4.2 Mean Field approaches 

Developing a consistent and well defined mean-field approximation of Eq.Q is 
not a trivial task [22|. It has been only recently that Birner et al. (completing some 
previous partial results [22]) have shown that [41 42 1 

A = sup (^T'l£)' <13) 

This result can be obtained by approximating the Laplacian in the Langevin equa- 
tion by its mean field analog, Y^Vi(4>(h t) ~ 4>( x , t)) = (4>) — 4>(x, t), writing 
the equivalent Stratonovich (see Appendix A) Fokker-Planck equation 

2 

dtPtf, t) = 84, (K + b4? - D({4>) - cf>)}p((f>, t)) + a — [<t>p(<f>, t)}) 

(14) 

and then solving it self-consistently in the stationary limit (see BT1 and |9l |22l 
l42l l. Let us underline the non-trivial nature of the mean-field result Eq. dT3l : for 
high noise-intensities it does not coincide with the zero-dimensional (3 (previously 
called Pi) exponent j3 = l/(p — 1), as usually occurs. Therefore, coupling to 
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the neighbors is a key ingredient to build up a sound mean-field description. At 
a field theoretical level (see next subsection) no systematic naive-power counting 
analysis or tree level approximation has succeeded so far in reproducing the non- 
trivial result of Eq.fHt (see \ 22\). Other critical exponents, can be easily computed 
as usual in mean-field theory. They are defined through r ~ A Uj - , £ ~ A"H , (4>) ~ 
t~ e , or r ~ £ z , with r (£) being the correlation time (length), and A = \a — a c \ 
the distance to the critical point. In the current mean field approximation: v±_ = 1, 
= 1/2, 9 = (3 = sup{l/(p - l),a 2 /(2D)}, and z = 2. 

4.3 Beyond Mean Field - Perturbative expansion 

In order to properly incorporate fluctuations and go beyond mean field, as usual, 
one has to resort to perturbative expansions combined with renormalization group 
(RG) techniques I24II431 . This can be done either directly at the Langevin-equation 
level, or alternatively working with the corresponding generating functional. Here, 
contrarily to [ 19] we choose to employ the second of these frameworks. The gen- 
erating functional associated with Eq.([l} 10 is l24l : 

Z[J,H] = V<j>Vt(]exp(-S[(t>,ip,J,H]) 
S = -Jd d xJ dt[^-(f) 2 ip 2 - + a<t> + b<pP - DV 2 (p] -J<f>- H<P] (15) 

where ip is a response field, J and H sources coupled to the order parameter field, 
(f)(x,t), and to the response field, ip(x,t), respectively 11 . Naive power counting 
shows that the upper critical dimension, below which mean field results cease to be 
valid, is d c = 2 12 . Following standard perturbative techniques, one can perform 
an e-expansion below d c . The noise amplitude is written in dimensionless form as 
a 2 = u[i e , where \i is an arbitrary renormalization scale (for example /i 2 = a). In 
momentum-frequency space the propagator is (—ito + a + k 2 ) (for simplicity 
we fix here D = 1) and the noise vertex is 




Here we use the Ito representation, as for it the Jacobian of the employed transformation and 
the theta function evaluated at can be seen to be equal to 1 and respectively, simplifying the 
algebra 1431 . The calculation would be almost identical in the Stratonovich representation, leading 
to essentially the same results, except for a shift in the critical point. 

"This functional can be easily worked out, by writing the probability of a given time-sequence 
realization of the noise as a product of Gaussian noise distributions at each time, generating in this 
way a path integral. Then variables have to be changed from r/(x, i) to 4>(x, t), and finally a dummy 
(ghost) variable ip is introduced by means of a Gaussian transformation, in order to avoid singularities 
typically induced by the MN. Afterwards the sources J and H are introduced by hand. The process 
is rather standard; a detailed description can be found in 1241 or l43l. 

12 This result is easily obtained by imposing the noise term coefficient to be dimensionless at d c . 
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The vertex (ilxj?) is like the previous one but with only one wavy line and p straight 
ones. It is easy to see that there are no perturbative corrections to the propagator, 
and the only possible diagrammatic corrections to the noise vertex have the follow- 
ing structure 



4 = - 2 



+ a 6 



a 



a 2 



I -a 2 



which can be written as o\ = jz~tt where I denotes the one loop diagram at zero 
external frequency and momenta, and /u 2 = a 



' f d d Kdu- % 5 — \ 7 = a^ (16) 



(27r)( d+1 ) J ioJ + k 2 + y 2 -iuj + k 2 + fi 

with a = 2/(8vr) d / 2 r(2 - d/2) 03^3] El- Th e pole 1/e captures the presence 
of an infrared singularity below d = 2. For d > 2, / is finite (i. e. it is infrared 
convergent) 13 therefore if a 2 is small enough a 2 1 < 1 and the mean-field theory 
predictions remain valid perturbatively. At a 2 = I a singularity appears in a 2 R ; 
above this noise value the perturbative expansion ceases to be meaningful; usually 
it is assumed that values of a 2 > I flow to a non-perturbative fixed point, con- 
trolling the long-wavelength large-time asymptotic behavior. For the corrections 
to vertex ipcf?, one has exactly the same structure, except for the rightmost noise 
vertex (see Eq. lfToT l) that is replaced by the other one, leading to 

b R = 6(1 + a 2 I + a 4 / 2 + ...) = YZ^2] = ha R/ a2 < 17 ) 

Therefore, renormalizing a 2 keeps also bn finite, and no extra renormalization is 
needed. As the degree of the non-linearity p does not introduce different diver- 
gences in the perturbative expansion, it should not affect the critical behavior (as 
long as p > 2). 

13 We assume the existence of a lattice cut-off ensuring the integrals to be ultraviolet convergent. 
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d<2 d=2 d>2 

DIMENSION 

Figure 2: Schematic representation of the RG flow of the dimensionless variable u 
as a function of dimensionality, as derived from Eq.Q. 



Proceeding with the renormalization program, we compute the Wilson (3— 
function defined as 

(3(u R ) = (id^aUR = -eu R - au 2 R (18) 

with a unique physical solution, ur = 0, stable for e < (i. e. above two dimen- 
sions) with associated mean field critical exponents. The renormalization-group 
flow-diagram is very alike the well-known one for KPZ P31I44 I. in particular, for 
d < 2 the flow runs towards infinity for any noise amplitude, while for d > 2 there 
is a separatrix between the trivial fixed-point basin of attraction and the runaway- 
trajectories, i. e. between the weak- and the strong-coupling fixed points. Let us 
emphasize both, the absence of any stable non-trivial perturbative fixed point, and 
the fact that the theory is super-renormalizable, in the sense that perturbative dia- 
grams have been computed and re-summed to all orders. An alternative derivation 
of these results can be found in [ 19 1. See also l44ll for a very similar calculation 
presented in a pedagogical way. 14 

So far, we have not proven that the strong coupling fixed points in MN (KPZ 
plus a wall) and KPZ coincide. 15 Observe, however, that right at the critical point 
the interface falls continuously to minus infinity, as far from the wall as wished. It 
is therefore very reasonable to assume, that asymptotically, the dynamical exponent 

14 It is illuminating to compare the generating functional Eq. <15> with that associated with the 
reaction A + A — > (kink-antikink annihilation) |46 47 32 48 1. It has the same structure as 1151 
but for the sign of the noise term. In terms of the Langevin equation the sign alteration tantamounts 
to the introduction of a purely imaginary noise 1321 12 |. This simple change has dramatic effects, 
in particular the Wilson /3-function is as Eq. <18t but with a positive sign for the quadratic term, 
generating a non-trivial perturbative fixed point below d = 2 [46 1, and a different non-trivial physics. 

15 It is well known that the KPZ equation should have a strong coupling fixed point controlling its 
non-trivial physics 1 441 1341 . 
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should be controlled by KPZ scaling (in the absence of the wall). In particular, it 
can be shown [20| that: (i) the dynamical exponent for MN1 coincides with that of 
KPZ, (ii) the correlation length exponent, v (that can be defined for MN but does 
not exist for KPZ, lacking of a critical point) can be related to the z exponent in 
KPZ by v = l/(2z — 2), and (iii) the order-parameter critical exponent (5 has to 
obey f3 > 1. The reader interested in the proof of these properties is deferred to 
[20 1 and to [ 30 1 for a more detailed explanation. 

4.3.1 Large N limit 

Exact RG calculations can be performed in the large N limit for the 0(N)-symmetric 
version of the MN equation defined as 

d t Mx,t) = -a<f> a - b\<j)\ 2 <j) a + DW 2 Mx,t) + a\<f>(x,t)\rj a (x,t), (19) 

where 4>(x,t) = {<fti,4>2, ■■■,4>n} is an N-component real vector field, and rj is 
a noise vector field. The cubic term is necessary in order to preserve the O(N) 
symmetry. The continuous symmetry of this model allows for two distinct types 
of active phases: one that preserves the symmetry and one which breaks it. The 
first one has (</>) = and Q = ([01 2 /N) ^ 0, while for the second both order 
parameters are non-zero. The main result of the exact solution presented in llT9l is 
that for d < 2, only the symmetric active phase exists, while for d > 2 both types 
of active phase can exist: for large values of a 2 the active phase is symmetric, 
while a broken-symmetry state emerges below a certain critical value of a 2 . There 
is also a multicritical point at which the two active and the absorbing phase join 
(see [19 1 for more details). At d = 2 a Kosterlitz-Thouless type of transition could 
appear 1 50 1. 

4.4 Microscopic models and Numerics 

Given the lack of analytical tools allowing for a determination of critical expo- 
nents in the strong coupling regime, one has to resort to numerical analysis. The 
first estimation of critical properties was obtained in ll20l by discretizing Eq.Q 
and performing a numerical integration in d = 1. A more complete study was 
presented afterwards in l22l in which the two- and three-dimensional cases were 
also analyzed. 

In what respect to discrete, microscopic realizations of this universality class, 
amenable to be simulated, let us remark that any of the many discrete models in 
the KPZ class reported in the literature 1331 l35l l34l is, in principle, susceptible to 
be modified by introducing a limiting wall, which induces a phase transition. If 
the non-linearity sign (that can be numerically determined by performing velocity 
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measurements of tilted interfaces [34]) is positive (negative) an upper (lower) wall 
is required to construct a model in the MN1 class. To the best of our knowledge, at 
least three different discrete models with these requirements have been studied in 
the literature H2"T1I5T1I52"1 . 

In finite size analysis of Eq.Q 16 one starts with homogeneous initial condi- 
tions, and first determines (</>) as a function of the system size, L, for different 
control parameter values. In a double-logarithmic plot the critical point corre- 
sponds to the straight line separating curves bending upward (pinned phase) from 
downward-bending ones (depinned phase). In order to have a well defined pinned 
phase for small systems one has to restrict the statistics to pinned runs (otherwise, 
the only true stationary state for finite systems is the depinned or absorbing one 
Il28ll29ll30l ). In systems with MN one has to define a criterion to decide when a 
configuration is depinned. Numerical results should not (and do not) depend on 
such a criterion. From the slope of log(^) versus log(i), one extracts 6 = 0/v\\ 
(see Fig.© to have a glance at how finite size scaling graphs look like). From the 
scaling of saturation values at criticality for different L one obtains (3/v±. Con- 
sidering a large system size it is possible to obtain (3 from the scaling of saturation 
values at various distances from the critical point. In what respect the dynamical 
exponent, z, it can be determined either from studies of the two-point correlation 
function at criticality [20] or from spreading experiments 12 1 1 1221 17 . Using spread- 
ing experiments at criticality the exponents 5 and rj, associated with the surviving 
probability and evolution of the total number of pinned sites respectively, can also 
be measured ll28ll29ll30ll3Tl. 

In Table 1 we list the values reported so far in the literature for d = 1. Even 
though there is some dispersion in the data, and a high precision determination 
of critical exponents is still lacking, the existence of universal values is firmly 
established. In d = 3 the two different universality classes predicted by the RG 
calculation (see Fig.©) can be easily recognized: for small a values very close to 
mean field expectations have been reported [22] while, for a larger a non trivial 
exponents have been found [22| (see Table 2). 

4.5 Realizations 

Among the physical realizations of this class, probably the most relevant one is a 
rather generic class of synchronization transitions in extended systems (see section 
lO Imi54ll52ll55ll56l. 

16 See 1 53 1 for a review of numerical methods designed to deal with stochastic equations. 
17 In spreading experiments 4>(x) = at time t = at every single sites, except for a seed-site. 
See (281129115(5). 
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Table 1: Critical exponents measured for the MN1 class in d = 1. RSOS stands 
for "Restricted solid on solid" model [52]. MN stands for simulations performed 
directly on the Langevin equation 
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Table 2: Critical exponents measured for the MN1 class in d = 3, for both the 
weak and the strong coupling regime 1221 . Simulations were performed directly 
in the Langevin equation with p = 2. The two different regimes are obtained by 
fixing the noise amplitude to a small and a large value respectively; see ll22l for 
further details. (3 m f is given by Eq. (fT3l . 



Another important application is in growing interfaces with both adsorption 
and evaporation of particles on a substrate. This may lead to pinning-depinning 
transitions or to non-equilibrium wetting transitions IBTI l57l l58l l59l (see I8.lt . 
It might also be related to the delocalization transition in an adsorption-reaction 
model reported in [ 60 1. 

In reaction-diffusion systems, many times it is possible to derive a Langevin 
equation in an exact way. In many cases, a multiplicative noise appears, but many 
times it is purely imaginary! (see the last section of this review refmix, and the 
references cited there), leading to a completely different phenomenology. 

5 The MN2 class 

As already anticipated keeping the KPZ non-linearity to be positive, the presence of 
a lower wall leads to a MN equation including an extra term with respect to the pre- 
viously studied case: Eq.©. It has a phenomenology substantially different from 
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that of MN1, Eq.Q. In this section we underline the analogies and differences 
between these two cases by going through the different levels of approximation 
discussed in the previous section. 

5.1 Zero dimensions 

The zero-dimensional case is obviously identical to the one for MN1. 

5.2 Mean Field approach 

Concerning the mean field approximation, the average value of the lower-wall ex- 
tra term, —2D(Vcl)(x,t)) 2 /(j)(x,t) generates an exponential non-integrable singu- 
larity, oc exp(+(0) 2 /</> 2 ), at the origin of the quasi-stationary potential, solution 
of the associated Fokker-Planck equation, rendering it non-integrable. Only the 
depinned (absorbing) phase exists within this approximation. A discussion of dif- 
ferent types of mean-field approaches, all of them leading to this same conclusion 
can be found in [62|. 

5.3 Beyond Mean Field - Perturbative expansion 

The additional term —2D(V4>) 2 /4> is non-analytical at <fi = 0, rendering un- 
feasible any type of perturbative expansion. Renormalization group results analo- 
gous to those presented for the MN1 class are therefore not available. The argu- 
ments leading to v = l/(2z — 2) obtained for MN1 can be extended to this case, 
but the one leading to (3 > 1 cannot. Indeed, as shown below, in this case v 1 as 
expected, and takes a value not far from 1/3. 

5.4 Microscopic models and Numerics 

In what respects numerical results, direct simulations of Eq.© are hopeless owing 
to the presence of the singular term. Contrarily, simulations performed in the KPZ 
representation are feasible, although they are also hard-to-deal-with as is geneti- 
cally the case when discretizing KPZ-type of equations [61 1. 

At the light of this barren perspective, one has to resort to simulations of mi- 
croscopic models. The first study was that of [21 1, where critical exponent close 
to those of the zero-dimensional case were reported, using a well known KPZ-like 
model I2T1 1351 . However, it has been recently shown that the results reported in 
II2TI correspond to effective exponent values computed in a transient regime , 

18 Apart from the numerical estimates, the rest of the discussion presented in 1211 about the differ- 
ences between the upper and the lower wall cases remains valid. 
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f3 u± f3/v±_ z 9 7] 



Hwa[62| 0.32(3) 0.97(5) 0.33(2) 1.55(5) 0.215(15) 

RSOS |j2) 0.325(5) « 1 0.33(2) « 1.51 0.21(1) 0.80(2) 



Table 3: Critical exponents for the MN2 class in d = 1, as measured in two differ- 
ent models: model 1 stands for the one described in 12T1 . and Model 2 is a RSOS 
one. See [62| for more details. 

and new accurate values have been measured (results are listed in Table 3 [62|). 
These studies are performed: (i) using a restricted solid-on-solid (RSOS) model 
l33l l34l l35l 152 1 known to belong, in the absence of a wall, to the KPZ class, and 
reaching system sizes as large as 2 20 , and (ii) the model considered in ETTl but 
with much larger statistics and performing finite size scaling analyses 19 . Results 
for both models coincide within errorbars as shown in Table 3. 

These values are clearly different from those reported for the upper wall, con- 
firming once again the dramatic (ever-surprising) difference between the two classes, 
as well as the existence of the robust MN2 class itself. 

The situation is much less satisfactory here than it is for the MN1 class. We 
have found problems both in mean field and in field-theoretical approaches, and 
even the numerical integration of the Langevin equation, in either the density of 
the interface representation is problematic l62l This calls for the development of 
new analytical tools aimed at tackling this problem and clarifying the physics of 
this (well identified and robust) universality class. 

5.5 Realizations 

The problem of probing the alignment of different DNA sequences is of outmost 
importance in molecular biology. One of the most widely used algorithms for such 
a problem is due to Waterman [36]. Hwa and collaborators, showed that such an 
algorithm can be mapped onto a KPZ problem with a lower wall (in particular, to 
its DPRM counterpart) [36 1. Locating the critical point of such a transition is what 
is effectively done by molecular biologists when tuning parameters ii. e. choosing 
a scoring scheme) to enhance or maximize fidelity in local alignment studies [36]. 

The MN2 class appear also in some pinning-depinning or non-equilibrium wet- 
ting transitions in which a KPZ interface (with positive non-linearity) grows from 
a substrate I5T1I2T1 . 

19 The discrete model considered in 1511 is extremely slow to simulate, as exhibits manifestly 
severe crossover effects. 
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It is known that one-variable MN processes under various circumstances (as, 
for example, the presence of reset events) lead to generic scale invariance [4(11 1^1 
l64l . The MN2 equation is the natural extension of this phenomenon to extended 
systems, and has been proposed, for instance, to describe the properties of city 
growth [63]|(3), where fluctuations in population number are multiplicative and, 
obviously, bounded from below. 



6 Changing the potential: MN1 with an attractive wall 

Up to now we have dealt with Langevin equations describing continuous transitions 
in systems with MN. In the same way as the Model A for Ising like transitions 
Eq.© has to be modified to the following Langevin equation 

d t <p(x, t) = -a<j> - b(p 3 - c4> 5 + DV 2 (j)(x, t) + ar](x, t). (20) 

with b < and c > 0, in order to describe discontinuous, first-order, transitions (in 
the absence of external magnetic fields) one might wonder what is the MN version 
of Eq. d20l and what type of physical problems could it represent. To explore this 
possibility one just needs to change the potential in Eq.Q to allow for discontinu- 
ous jumps: 

d t 4>{x, t) = + DV 2 <j)(x, t) + a<t>{x, t)rj(x, t) 

o<p 

V(0 = ^ 2 + ^ 3 + -/ (21) 

with b < and c > 0. If (j) — (j) symmetry is to be respected, the simplest 
potential becomes acj) 2 + 60 4 + cc/) 6 (with no essential difference in the underlying 
physics). Transforming Eq. d2"Tt to the interface language: 

d t h(x, t) = + £>V 2 /i + D(Vh) 2 + arj(x, t) 

V(h) = ah + be h + c/2e 2h . (22) 

For b < this is a KPZ equation in the presence of an attractive upper wall, i. 
e. in order to depin a site close to the wall, a certain potential barrier has to be 
overcome. Indeed, the potential V(h) has a local minimum nearby the wall (see 
Fig-®- In analogy with the previous sections, this KPZ equation with a positive 
non-linearity and an upper wall is fully equivalent to a negative non-linearity equa- 
tion with an attracting lower wall. Eq. d2TT > is therefore the b < counterpart of 
the MN1 class. This equation has been studied in different contexts in the litera- 
ture. It was first proposed by Miiller et al. [ 66 1 to study first order non-equilibrium 
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Figure 3: Effective potential in terms of h for different parameter values. For b < 
it has a minimum nearby the wall. 



transitions in extended systems. Afterwards, a very similar equation was proposed 
by Zimmerman et al. ll67l as a simplified stochastic model of spatio-temporal in- 
tennittency, and by Giada and Marsili f42l to study the effect of Morse potentials 
in KPZ dynamics 20 . Also, in an interesting paper, Hinrichsen et al. ll68l studied 
what they called first order non-equilibrium wetting: a discrete random deposition 
model, known to be described by KPZ, in the presence of both evaporation (i. e. a 
drift) and a limiting attractive substrate (i. e. the motion of the KPZ interface is not 
just limited by the wall, but the wall is "sticky" in the sense that sites attached to 
the wall need to overcome an extra potential-barrier in order to become depinned) 
was analyzed. It can be easily argued [68 ] that this discrete model is represented at 
a continuous level by Eq. d2Tl l. 

Following the same steps as in preceeding systems we now analyze Eq. dTiTi 
employing different methods. 

6.1 Zero dimensions 

In zero dimensions, the problem is exactly solvable. The stationary probability dis- 

tribution being (see appendix B) P s t\4>) = ^ 1+2a/tT ^ e a \ > . At the transition 
the order parameter jumps discontinuously from to a non- vanishing value. 

6.2 Mean Field approach 

At a mean field level, we can perform a calculation analogous to that presented 
for Eq.©. The numerical solution of the corresponding self-consistency equation 

20 Observe that V(h) has the form of a Morse potential. 
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Figure 4: Phase diagram derived from the numerical integration of the mean field 
approximation of Eq. d2"TY In the coexistence regions both, the pinned and the 
depinned solutions are stable, and their edges (marked by dashed lines) correspond 
to discontinuous transitions, i. e. one of the phases loses its stability abruptly. The 
transition corresponding to the upper one, a c (b) (signaling the limit of stability of 
the pinned phase), may become continuous once fluctuations are considered, as 
discussed in the text. 



generates the phase diagram shown in Figure ©. We describe now the main fea- 
tures (for a more detailed explanation see l66l l42l BTl l58lD . For positive values 
of b the transition occurs at a = 0, it is continuous and very similar to the one in 
Eq.Q: it is in the MN1 class 21 . Instead, for b < the transition bifurcates: there 
is a discontinuous transition at a = 0, and a second transition (also discontinuous 
at mean field level) appears at a positive value of a (see Fig.©), a c (b). At this 
value, the pinned phase, (<p) ^ 0, loses its stability and falls in an abrupt way to 
the depinned phase = 0. The depinned phase loses its stability at a smaller value 
a* = 0. Therefore, in the interval [0, a c (b)] both solutions are stable: we are in 
the presence of phase coexistence. The existence of broad coexistence regions is a 
trait specific of non-equilibrium phase order transitions, and have been reported to 
appear, among others, in the Toom model [69]. Observe that both of the limiting 
lines of the coexistence region are discontinuous transitions at this level. 

6.3 Beyond Mean Field - Perturbative expansion 

The theory can be treated perturbatively in a very similar way to Eq.Q; indeed 
the Feynman diagrams have the same topology in both cases. The new constant 
c renormalizes as cr = c(l + a 2 1 + a 3 1 2 + ...) = 1 _ c a '2 I and a RG flow anal- 
ogous to that for Eq.Q is obtained. The only difference is that now as long as 

21 As in MN1, in this case, for a > (a < 0) the tail of the potential would attract (repel) distant 
interfaces to (away from) the wall; the transition from one case to the other occurs at a = 0. 
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the transition is discontinuous (not scale invariant) the trajectories for large enough 
noise amplitudes (above the separatrix in Fig.©) are not expected to flow to a 
strong coupling fixed point, but to be runaway trajectories. In any case, there is no 
non-trivial stable fixed-point and, consequently, standard perturbative expansions 
cannot say anything about the strong noise regime. 

6.4 Microscopic models and Numerics 

Eg. dm has been studied numerically in different contexts [67 42 1 and recently in 
l57l l58l . Its main properties (described mainly in the interface language) are as 
follows. The mean-field phase-diagram qualitative-structure survives the inclusion 
of fluctuations, namely: 

6.4.1 b > 

For b > a continuous phase transition is observed. This transition has been 
extensively verified to belong to the MN1 class l57ll58l . 

6.4.2 b < 

For b < simulations performed taking D = 1 confirmed the existence of a broad 
coexistence region. Within this region an initially detached interface falls towards 
minus infinity with a finite average velocity (i. e. absorbing remains absorbing) 
while an initially pinned region remains pinned. This is strictly true only in the 
thermodynamic limit. For finite sizes, the only stable state is the depinned (absorb- 
ing) one l28ll30ll . The characteristic time required for an initially pinned interface 
to detach goes to infinity for large system sizes within the coexistence region, con- 
firming the thermodynamic stability of the pinned phase. Deep into the pinned 
(active) phase this characteristic time grows exponentially, while close to a c {b) 
there is a tiny region where such a time seems to grow as a power law 22 . Within 
this regime, spatio-temporal coexistence of the two phases is actually observed in 
computer simulations 23 . In Fig. © we see a space-time plot of the activity in 
this regime. Pinned and depinned regions coexist implying that there are stable 
mechanisms for both: generating fluctuations and eliminating islands of the minor- 
ity phase. In this way a spatio-temporal intermittency (STI) regime [ 70 1 is shown 

22 In the absence of a theoretical understanding, it could well be the case, that this power law is 
just a transient, becoming asymptotically an exponential. Indeed this seems to be the most likely 
scenario. 

23 Noteworthy that phase coexistence (between a pinned and a depinned phase) is one thing, and 
spatio-temporal coexistence of detached and non-detached patches is a different one. This last is 
usually identified with spatio-temporal intermittency. 
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Figure 5: Spatio-temporal evolution of Eq. d2Tl in the spatio-temporal intermittency 
region, close to a c (b). In the horizontal axis we plot the time, while the space- 
position is plot in the vertical one. Red (light) for pinned sites, i. e. large (small) 
(j), and blue (dark) for depinned ones. 



to exists l67l BTl . Analogous mechanisms, required to sustain two-phase spatio- 
temporal structures, have been described in detail in a microscopic realization of 
this class [68 71 1. Performing a cut for a given time, and plotting h = — log(</>) 
as a function of the position, we obtain a profile as the one shown in Fig. (0. The 
depinned regions tend to escape from the wall, but being pinned at their edges, 
they are constrained to form very characteristic triangular patterns. The charac- 
teristic slope, s of the triangles can be determined by observing that A((V/i) 2 ) has 
to compensate the effect of the constant drift a, and therefore As 2 = a. As A can 
be measured from the average speed of tilted interfaces 1341 1331 . s can be inferred. 
Indeed, the value obtained in this way is in very good agreement with the measured 
average slope as reported in [58][72- Summing up, inside the coexistence region 
we have two different sub-phases. In both of them the depinned phase is stable, 
but: 

1 . One has small fluctuations of the pinned phase, and exponential increase of 
detaching times Two homogeneous phases, pinned and depinned, coexist. 

2. The other has larger fluctuations in the pinned phase, allowing for finite 
patches to depin and form typical triangular structures. In this regime spatio- 
temporal intermittency is observed. The detaching times in this regime seem 
to grow algebraically with L although this could be a transient effect. In this 
regime a STI-pinned phase coexists with the homogeneous depinned one. 

In higher dimensional systems the role of the triangles is played by pyramidal 
structures (see Fig. 0and [58].) 
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Figure 6: Values of /i = — ln(^>) in the previous figure, at a given fixed time-slice 
(roughly speaking at 3/4 of the maximum time; the large triangle, corresponds to 
the up-right absorbing region in the previous picture). 




Figure 7: Two views of a pyramidal structure appearing in two-dimensional simu- 
lations, in the phase coexistence region, of the MN noise equation in the presence 
of an attractive wall. 
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6.4.3 A DP transition 



Surprisingly enough, numerical studies have recently shown that there is a region 
in parameter space for b < 0, not predicted by the mean field approach, in which 
the transition at a c (b) becomes a continuous one. We will illustrate in what follows 
that the presence of fluctuations may induce the mean-field first-order transition 
at which the pinned phase looses its stability to become second-order, governed 
by DP exponents. This is somehow surprising, as the broad coexistence region 
survives, but its upper border edge corresponds to a continuous transition. 24 

In [ 18 1 it has been shown that, at least, for certain parameter values, satisfying 
the following two requirements: i) to generate at mean field level a weak (small 
jump) discontinuous transition (in such a way that fluctuations can alter more easily 
its degree), and ii) to have a relatively small coupling constant D, so the mechanism 
at the basis of the instability of the pinned phase is weakened, 25 a continuous DP 
transition is generated. The detailed phase diagram has still to be worked out and 
the conditions under which a continuous transition is generated need to be further 
clarified. 

What is already firmly established from numerical simulations in fl"8l is that 
such a continuous-transition regime exists and that it is controlled by DP expo- 
nents. In particular, in Figure [8] we show the order parameter (</>) time evolution 
for different system sizes. Both the slope at criticality and the scaling of satura- 
tion values for finite sizes, show unambiguously DP values (see figure caption). 
In order to fully verify the presence of DP scaling, different tests have been per- 
formed in [18]. In particular, the moment ratios have the same (universal) values 
as those measured for the DP-class [74]. It was also verified that the number of 
sites pinned by the potential scale with the same (DP) exponents. Finally, by com- 
puting numerically the fluctuations from different initial states, it was checked that 
they are proportional to the square-root of pinned (active) sites, as occurs in the 
Reggeon-field theory describing DP I28l l30l . 

To further confirm the emergence of DP scaling with b < 0, and in order to 
confirm the generality of the described mechanisms, a well-known discrete inter- 
face model known to be in the KPZ class (i. e. the one described and studied in 
Eli ) has been studied by including an attractive wall. The introduction of attrac- 
tiveness changes in a very abrupt way the critical behavior: while in the absence 

24 A similar situation (second order transitions with first-order traits) named "one-way hysteresis", 
has been recently reported to occur in some different contexts. See R. Maimon and J. M. Schwarz, 
condmat/0301495 and references therein. 

25 Detached triangles pulling with a strength which depends on D are responsible for the destabi- 
lization of the pinned phase 1 18 68 1 (See |71 1 for a nice explanation of the triangles dynamics and 
stability. 
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Figure 8: Log-log plot of the order parameter time evolution (averaged over 
surviving, pinned, runs) at the critical point, for system sizes (top to bottom) 
L = 25,50,100,200, 400, and 1000. Both the main slope, 9 = 0.159(5), and 
the scaling exponent of saturation values, (3/v± = 0.245(15) reveal DP-scaling 
behavior (0 = 0.15947... and f3/u ± = 0.25208(5) |73|). Inset: Log-log plot 
of the average order parameter as a function of \a c — a\. The full line is a fit 
(<fi) ~ \a c — a\P, with an exponent J3 = 0.28(1), also compatible with the DP value 
(3 = 0.27649(4). 



of attraction (purely bounding wall) the transition is described by the MN1 class, 
as soon as attractiveness is switched on (and therefore b changed from positive to 
negative) the transition crosses over to DP. Combining Monte Carlo simulations 
with finite size scaling, many decades of clean DP scaling can be observed. 

According to a recent preprint by Droz and Lipowski [75], given that: i) the 
absorbing (depinned) state is compatible with multiple possible realizations, and ii) 
at the light of numerical simulations, this transition is in the DP class with infinitely 
many absorbing states [76]. This means in particular that the spreading exponents 
may be not universal 1761 . 

Before ending this section, let us remark that it could well be the case, that this 
transition is controlled asymptotically, in d = 1, by a DP fixed point for any param- 
eter values. Indeed, some heuristic arguments have been proposed by Hinrichsen 
pointing out to the impossibility to have true first-order transitions between an ab- 
sorbing and a fluctuating phase in one-dimensional systems l49l . If this scenario 
is true, then the apparent discontinuous transition at the upper border line of the 
coexistence region would be just a transient, crossing over to DP. In this scenario, 
the described STI phenomenology would corresponds to DP-dynamics. In higher 
dimensions, first order transitions can appear without troubles. 
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6.5 Open problems 



An open problem is how to deduce analytically (i. e. from a RG perspective) that 
Eq. d2Tl can generate DP scaling for some parameter values. 

An intuitive, heuristic justification of this fact is as follows: the scaling, domi- 
nated in general by interface fluctuations is controlled here (in the DP regime) by 
fluctuations of the number of pinned sites. The attractive wall, by generating a min- 
imum in the potential, provides a physical ground to distinguish between pinned 
and not-pinned sites. We can assign them, I's or O's respectively, and define in 
such a way an effective dynamics for these booleans variables. This has to be very 
similar to DP. For instance, depinned (absorbing) regions will not become pinned 
(active) as, once away from the potential well, the interface is pushed exponen- 
tially towards infinity, so 1 can go to 0, but the reverse is not true. However, from 
a more rigorous point of view, it has not been possible so far to derive DP-scaling 
from Eq. d2TT> using RG techniques. Indeed, as anticipated previously perturbative 
studies can say nothing about the strong-noise regime allegedly controlled by DP- 
scaling 26 . 

Another problem is related to the possibility raised by Hinrichsen [49| that 
the first-order transition at which the pinned phase becomes unstable, is only ap- 
parently discontinuous. It might well be the case that the line of transition a c (b) 
is continuous all the way, with asymptotic DP-scaling. This needs to be further 
clarified in more extensive computer simulations. 

6.6 Realizations 

The main two realizations known so far of this type of MN equation in the presence 
of an attractive wall are non-equilibrium wetting and synchronization in extended 
systems. These two realizations are discussed in detail in forthcoming sections 

It could also be related to the Pair Contact process with diffusion (one of the 
most intensely studied and debated model with absorbing states [77 1) for which a 
similar equation with MN has been derived If78ll79l . as pointed out very recently 
by Hinrichsen [71 1. This could justify the conclusion that this controversial model 
exhibits asymptotically DP scaling. 

26 A very similar problem emerges in the context of fluctuating (KPZ) interfaces in the presence 
of quenched disorder. Also in that case, DP-scaling emerges in a way unaccessible to standard 
perturbative methods (see Appendix C). 



26 





MNl— transition 










•a 


Depinned 


« 




a 

a . 


i Tricritical Phase 


s 


\ Point 




\^ DP or First order. 


First 




Order 






Coexistence \ 




Region 





MN2 — transition 


iase 








Depinned 






Pinnei 




Phase 

► Tricritical Point 

First 
Order 



Figure 9: Schematic representation of the phase diagram for upper (lower) and 
lower (upper) walls with a positive (negative) non-linearity. See the caption of 
Table 4 
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Table 4: Universality classes for the KPZ equation with different signs of A, in 
the presence of Upper or Lower walls, with (YES) or without (NO) attractiveness. 
MNl stands for the universality class described by Eq.Q; MN2 stands for Eq.©; 
1 st represents discontinuous transitions, DP transitions in the directed percolation 
class, and BPC for a broad region of phase coexistence. In the case of a BPC 
region there are two phase transitions, one (where the depinned phase looses its 
stability) is always discontinuous (not reported in the table) and another (at which 
the pinned phase becomes unstable) which may be either first order (in dimensions 
larger than two, and probably in an apparent, transient way in d = 1) or DP-like, 
as shown in the table. Also, in this case, if the value of |6j is small enough, the 
renormalized value of b might be positive leading to a standard continuous MNl 
transition, without BPC. In the case of MN2 with attractiveness no BPC appear; 
the transition is second order, in the MN2 class, except for large negatives values 
of b where it becomes discontinuous. 
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7 The remaining possibility: MN2 with an attractive wall 



The equation described in the previous section is, as explained, the b < (attractive 
wall) counterpart of the MN1 class. To exhaust all the possibilities, we just need to 
explore the phenomenology of attractive walls (6 < 0) added to the MN2 class, i. e. 
lower attractive walls with positive non-linearities or equivalently upper attractive 
walls with a negative KPZ non-linear coefficient 27 . 
In the interface language we have 

d t h{x, t) = + DV 2 h - D(Vh) 2 + ar](x, t) 

V(h) = ah + b exp /i + c/2 exp(2/i) (23) 

with b < 0. This can be mapped into 

(Vd>(x t)) 2 

dMx, t) = -a<p - b(f> 2 - c<f) 3 + DV 2 <f>(x, t) - 2D y Y ; ' " + a<p(x, t)rj(x, t). 

<p{x,t) 

(24) 

with a negative value of b. Having arrived almost exhausted to this fourth possibil- 
ity we will not split this section into the usual subsections of the preceeding ones. 
We just stress that b < does not have effects as dramatic as it had in the previ- 
ous case. Here the transition, which at mean field level is discontinuous, becomes 
continuous for relatively small values of |6| (and belongs to the same universal- 
ity class as the b > case i. e. to the MN2 class.), while for sufficiently large 
ones, it remains first order without a broad coexistence region 28 . In this way, the 
presence of fluctuations may drive the transition from discontinuous (at mean field 
level) to continuous in some parameter range. In other words the tricritical point is 
shifted from b = at mean field, to a negative value when fluctuations are taken 
into account (this effect is much stronger here than it is in the previous case). This 
result has been verified in numerical simulations both using the model in [21 1 and 
a RSOS model in the presence of an attractive wall ll62l . Let us stress again that 
even if the transition is discontinuous, first-order, in the strong attraction regime, it 
does not involve a broad coexistence region nor a region of STI as occurred in the 

27 We have to admit that this looks like a tongue-twister and we humbly suggest the long suffering 
reader to draw a diagram to clarify the list of possibilities. This diagram should include 2 3 situations: 
upper or lower wall?, positive or negative non-linearity?, attractive or non attractive wall?; and only 
4 different physical behaviors, described respectively by this and the 3 preceeding sections. Table 4 
and Fig. j5J should be of some help. 

28 Observe that the first-order transition here does not require the presence of coexistence between 
active and absorbing phases: it just corresponds to the point where the free interface has zero velocity, 
and therefore does nor contradict Hinrichsen's claim 1 49 1. 
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previously discussed case. The transition is located for any value of b at the point 
where the depinned phase (free interface) looses its stability 

One way to rationalize this result (and its being at odd with the previous section 
one) is as follows. In the MN1 class in which there is an effective repulsion from 
the wall as discussed at the beginning of the paper, the introduction of attractiveness 
has a profound effect. In particular, two opposite physical mechanisms compete: 
one repelling sites from the wall and the other attracting them to it. This is at the 
origin of the existence of a broad coexistence region. Contrarily, for the MN2 class 
such an effective repulsion does not exist and, therefore, introducing attractiveness 
is a much milder perturbation in this case. Indeed, it does not alter nor the order 
nor the exponents of the transition that remains in the MN2 class for sufficiently 
small values of \b\, and it induces a standard discontinuous transition (no broad 
coexistence region) for large enough (very strong) attractiveness. 

A similar phase diagram, including first and second order transitions as well as 
a tricritical point separating them, has been reported by Hinrichsen et al. in their 
discrete non-equilibrium wetting model, for parameter values corresponding to the 
case under consideration. 



8 Applications: Non-equilibrium wetting and synchroniza- 
tion 

In this section we briefly review two of the main fields in which the formalism and 
results of the previous sections can be straightforwardly applied: non-equilibrium 
wetting and synchronization in spatially extended systems. We discuss how MN 
equations as those described in the preceeding sections come out within these con- 
texts. 

8.1 Non-equilibrium Wetting 

When a phase (a) is in contact with a substrate, wetting occurs if a macroscopic 
layer of a coexisting phase (/3) is adsorbed at the substrate. The wetting transition 
is characterized by the divergence of the /3-layer thickness; the a/5 interface goes 
arbitrarily far from the substrate. Equilibrium wetting is a problem of great experi- 
mental as well as theoretical importance. It has been extensively investigated using 
interface displacement models, in which h(x,t) represents the interface distance 
from the substrate 1 80 1 . The minimal dynamical (equilibrium) wetting model is 
given by the following Edwards-Wilkinson (EW) [81] Langevin equation [ 82 1 

dV 

3t/i(x, t) = DV 2 h -^r + <77/(x, t). (25) 
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This was introduced by Lipowsky 1821 . and describes the relaxation of h towards 
its equilibrium distribution, specified by V(h). A physically motivated choice for 
the interaction potential has the following form (Morse potential) ll80l 



where b(T) vanishes at the wetting temperature and c > 0. At sufficiently low tem- 
peratures, b < 0, the potential V binds the a(3 interface; the equilibrium thickness 
of the wetting layer (h) is finite. As the temperature is raised, the potential becomes 
less attractive and at some point it no longer binds the interface, (h) diverges. 

In order to study complete wetting (i. e. wetting occurring in the limit in which 
the chemical potential difference between the two phases, fj,, goes to zero for T 
larger than the wetting-temperature) a linear term fih has to be added to V(h) [ 80 1. 
\i plays the role of an external force acting on the interface. 

The dynamic model, Eq. d25l is readily generalized to non-equilibrium inter- 
facial processes, e.g., crystal growth, atomic beam epitaxy, etc., where thermal 
equilibrium does not apply. This has been done in detail in I571I581B91 . 29 The 
conclusion is that the most basic effective non-equilibrium interfacial model con- 
sists of a KPZ equation (instead of Eq. d25l ) in the presence of a substrate, i.e.: the 
MN equation. The origin and nature of the KPZ non-linearity within this context 
has been discussed in I57ll58ll59l . 

A detailed discussion about non-equilibrium critical wetting, complete wetting, 
the difference between pinning-depinning and wetting-dewetting transitions, and 
other related issues can be found in [ 58 59 1 (see also [ 83 ]). 



8.2 Synchronization 

Mutual synchronization of oscillators is one of the most intriguing and fascinating 
phenomena appearing in complex systems [84 85 1. Some examples appear in 
chemical reactions [86], neuronal networks [87|, flashing fireflies [88|, Josephson 
junctions [89|, and semiconductor lasers [90] to mention but a few. See [84] for an 
elementary introduction to this growing field and [85] for an excellent review. 

Coupled map lattices (CML) 1911192 ;! have attracted much attention as a paradigm 
in the study of synchronization in spatially extended systems under mathematical 
grounds. When two replicas of a same CML are locally coupled [55j|93l, or when 
they are coupled to a sufficiently large common external random noise even if they 
are not coupled to each other [94]|95]|54), they can achieve mutual synchronization. 

In all the aforementioned examples, there is a transition from a chaotic or un- 
synchronized phase, in which the differences between both replicas persist, to a 

29 See also I51II68I where the same problem was first studied for discrete models. 



V(h) = b(T)e- h + ce 



-2h 



(26) 
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synchronized one in which memory of the initial differences is erased and replicas 
synchronize with certainty. The analysis of universal properties of synchronization 
transitions (ST) has been the subject of many recent studies [55]|56l|54]EllI21- m 
particular, as mentioned previously, in a seminal paper, Pikovsky and Kurths de- 
vised a method to map the synchronization error field (i. e. the difference between 
two replicas) into the MN equation ^ On the other hand, it was proposed that for 
discrete systems the ST should be in the DP class I55l l56l . 

Based on recent observations ll54l l52l l95l a global picture, that can be synthe- 
sized as follows, has emerged: there are two types of ST depending on whether 
the largest (transverse) Lyapunov exponent, A, is zero or negative right at the tran- 
sition l54ll52l . In the first case, the transition is MN like, and it is controlled by 
linear stability effects (i. e. A changes sign at the transition). In the second group 
the transition is either DP or discontinuous depending on microscopic details (and 
the transition is located where the propagation velocity of non-linear perturbations 
changes sign, while A is negative). 

This situation is in complete analogy with what reported for Eq. d2TT i. Indeed, as 
shown in [18|, a straightforward generalization of the derivation by Pikovsky and 
Kurths leads to Eq. d2TT> instead of Eq.Q in the cases when there is some singularity 
or discontinuity in the local map. 30 In particular, as discussed in [ 18 1 if there is 
a typical re-injection rule or another mechanism leading to a typical value of the 
local synchronization error field, then the corresponding potential in the Langevin 
equation needs to have a local minimum, and therefore, the structure of Eq. d2TT> 
rather than Eq.Q- 

In conclusion, Eq. d2Tl is the minimal Langevin equation describing ST in a 
general way. All the so-far reported universality classes in one-component syn- 
chronizing system are describable by this equation, and consequently belong to 
one of the three classes reported for it: MN1, DP or are discontinuous. See lTl~8l 
for further details and analysis. 

9 Miscellaneous comments 

We present some miscellaneous remarks on issues related to the ones discussed 
along this paper; they might be of some interest in order to get a global picture of 
the state of the art. See also the recent review 1971 . 

• Negative Diffusion 

30 In a recent paper 1 96 1 it has been shown how DP can emerge in synchronization transitions, even 
for some continuous maps. 
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By considering negative values of the diffusion constant, anti-ferromagnetic 
ordering can be studied. In particular, Birner etal. [41 1 have recently studied 
a MN equation with p = 3. They consider both positive as well as negative 
values of the order parameter field and report on an anti-ferromagnetic MN 
phase transition (in the MN1 class). 

Additive Noise 

The presence of additive noise in an otherwise MN changes the universality 
of the transition (9] El [T6l : it is a relevant perturbation from the RG point 
of view. This is easily checked as the additive noise is more relevant from 
naive power-counting analysis than the multiplicative one. The phenomenon 
of reentrance may appear when both multiplicative and additive noises act 
together: the system can become more ordered by increasing the MN ampli- 
tude, but eventually the additive part will take over, disordering the system 
lfT31IT6l . For example, it has been shown that annealed Ising models, known 
to have a reentrant transitions can be described by Langevin equations with 
competing additive and multiplicative noises [98 1. 

Colored Noise 

One can wonder what is the effect of colored noise in MN type of equations. 
From the RG point of view deviations from Gaussianity in the noise are 
irrelevant, but they might have some effects on the phase diagram (see for 
instance IH [lOOl). In (T01| the effect of spatially oscillating MN has been 
studied. 

Higher order Noises 

Having studied Langevin equations with minimal potentials and noises am- 
plitudes proportional to 4>°, 1 / 2 , and cj) 1 , one might wonder if higher powers, 
as for instance 3/2 or 2 have any physical significance. If one construct us- 
ing standard Fock space techniques [48 1 an effective Langevin equation (or 
equivalently a generating functional) for reaction diffusion process occurring 
in n-plets (triplets, quadruplets, etc), it is easy to see that they generate noise 
amplitudes proportional to 4> n / 2 If all possible reactions in a given system 
require the presence of triplets, it is reasonable to expect a 3/2 type of noise 
(79] [102 1. For a general n — plet process a noise amplitude proportional 
to n / 2 is expected to be the leading contribution. This could generate new 
universality classes not fully explored yet. See also [103 1 where other noise 
terms are studied. 

Long range interaction 
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If one considers long-range (for example, van der Waals) interactions be- 
tween the substrate and the interface in the wetting problem the potential has 
the general form fT04ll80l 

V(h) = b(T)h- m + ch~ n , n>m>0. (27) 

A similar interaction has been recently used to describe a class of synchro- 
nization problems for which a new universality class has been claimed to 
appear [ 105 1. We are presently working in this problem. 

• Quenched disorder 

The presence of quenched disorder is known to generate new non-trivial phe- 
nomenology in both Ising and DP -like systems. It is also expected to be a 
relevant perturbation at the different MN universality classes reported here. 
Its analysis and relation with physical processes, as localization, is currently 
under investigation. 

• Multicritical Points 

At the points where different critical lines meet, a different, multicritical 
behavior is expected to show up. These multicritical points and their corre- 
sponding Langevin representations (most likely including MN) are still to be 
studied in depth [52 105 1. A multicritical point like these should be related 
to critical non-equilibrium wetting as described in l58l (see also [83 1). 

• Complex Multiplicative Noise 

It is well known that certain (annihilation) processes can be represented by a 
Langevin equation with MN, but a complex amplitude l46l l47ll48l l2l l32M78l . 
The degree of validity of such Langevin equations with complex noise still 
lacks of a satisfactory analysis; see [32, 78] for some discussions of this and 
related topics. 

10 Conclusions 

We have reviewed different aspects of non-equilibrium phase transitions exhibited 
by systems amenable to be represented by a Langevin equation with MN. These 
systems exhibit a phase transition from an absorbing phase to an active one, as 
for example pinning-depinning transition, non-equilibrium wetting transitions and 
synchronization transitions. A rich phenomenology can be further explored by ex- 
ploiting the Cole-Hopf transformation that maps MN systems into non-equilibrium 
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(KPZ) interfaces in the presence of an extra limiting wall. Depending on the po- 
sition and nature of the wall we have identified and characterized up to 4 different 
types of critical behavior that are synthesized in Table 4. 

It is our hope that this review will foster new developments in the fields of non- 
equilibrium phase transitions and systems with MN, and stimulate the application 
of these ideas to other open problems in physics and in other disciplines. 

Appendix A: Ito versus Stratonovich 

Let us consider a MN equation in the Stratonovich representation. Then it is 
completely analogous to an Ito-interpreted Langevin equation in which the linear 
coefficient a has to be replaced by a + a 2 / 2 UE). Therefore, for this type of noise 
the only difference between Ito and Stratonovich tantamounts to a shift in the con- 
trol parameter. Along this paper we use both Ito and Stratonovich interpretations, 
as the results obtained for any of them can be straightforwardly modified to apply 
to the other. An exception are noise induced transitions as those reported in [9j, the 
difference between Ito-Stratonovich is essential in this case, as the noise amplitude 
plays a key role in the linear control parameter only if the equation is interpreted 
a la Stratonovich fl6l . See [ 106] where a new type of noise induced transition is 
reported to occur for both Ito and Stratonovich interpretations. 

Appendix B: The 0-dimensional case. Solution of the sta- 
tionary Fokker-Planck equation. Comparison with RFT. 

Let us consider Eq.Q in the zero-dimensional limit 



The associated Fokker-Planck equation in the Stratonovich representation is 



whose stationary solution (trivially obtained imposing detailed balance) is given 
by Eq. lfTTIi . As long as a > the probability becomes non-normalizable, and the 
only stationary state is = 0. Therefore the transition is located at a c = 0. In the 
Ito interpretation it is shifted to a = a 2 /2. Observe that in the absorbing phase the 
singularity at the origin is non-integrable and changes continuously its degree but, 
above the critical point, the singularity disappears. Instead, for the RFT noise (i. e. 
replacing eft by 1 / 2 in the noise term), following the same steps one obtains, 



d t cp{t) 



a<p — b(fr p + a(j)r] 



(28) 



dtP(<t>,t) = d^[{acp + b^)P{4>)] + —d^cpd^P^)]] 



(29) 




(30) 
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Here the singularity exists in both the active and the absorbing phase, therefore 
P{4>) reduces to a delta function at the origin. This reflects the fact that for RFT- 
like finite systems the only true stationary state is the absorbing one. Different 
plots, comparisons and discussion of the (pseudo)-stationary potentials for RFT 
and MN can be found in [ 32 1. 

Appendix C: The problem of the quenched KPZ 

Let us consider the problem of a KPZ interface, with a positive non-linearity 

in the presence of quenched disorder (QKPZ), h), where £ is a Gaussian dis- 
tribution. It was proposed in [107] that the scaling properties of this problem at 
criticality should be controlled by a DP type of fixed point in d = 1 3 . The logic 
behind this prediction is that the interface is pinned by a path of impurities such 
that it percolates from one side to the other, blocking the whole interface front. As 
the path is directed, DP is expected to control its scaling. Observe that, owing to 
topological reasons the same reasoning does not apply to two dimensional systems. 
In d = 2 a directed surface is required to block the interface advance. A scaling 
theory of directed surfaces was proposed in [ 108 1 to account for this. 

Despite of the numerous computational confirmations of the QKPZ-DP con- 
nection, no analytical study has succeeded so far in establishing such a link in a 
more rigorous way. The difficulties stemming from the fact that it seems that the 
order parameter has to be changed: it is the interface velocity for KPZ, while it 
should be related to the density of pinning sites in the QKPZ case. Establishing 
under firm grounds the equivalence between QKPZ and DP remains an open chal- 
lenging problem. 

In the case of a negative KPZ non-linearity, simulations show either triangu- 
lar patterns resembling the ones discussed previously 32 and a first order pinning- 
depinning transition has been reported to occur f72l or, alternatively, critical prop- 
erties controlled by DP [ 109 1 depending upon model details or parameter values. 
Let us emphasize the deep analogies between this situation and the one described 
in Section |6] 
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